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ABSTRACT 

We study the properties and stability of nuclei in the inner crust of neutron 
stars in the presence of strong magnetic fields ~ 10^^ G. Nuclei coexist with a 
neutron gas and reside in a uniform gas of electrons in the inner crust. This prob- 
lem is investigated within the Thomas-Fermi model. We extract the properties of 
nuclei based on the subtraction procedure of Bonche, Levit and Vautherin. The 
phase space modification of electrons due to Landau quantisation in the presence 
of strong magnetic fields leads to the enhancement of electron as well as proton 
fractions at lower densities ~ 0.001 fm~'^. We find the equilibrium nucleus at 
each average baryon density by minimising the free energy and show that ,m the 
presence of strong magnetic fields, it is lower than that in the field free case. The 
size of the spherical cell that encloses a nucleus along with the neutron and elec- 
tron gases becomes smaller in strong magnetic fields compared with the zero field 
case. Nuclei with larger mass and atomic numbers are obtained in the presence 
of strong magnetic fields as compared with cases of zero field. 

Subject headings: stars:neutron - inner crust - nuclei - magnetic fields 



1. Introduction 

Strong surface magnetic fields ~ 10^^ G are found to exist in pulsars. Even stronger sur- 
face magnetic fields > 10^^ G were predicted by observations on soft gamma ray repeaters and 
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anomalous x-ray pulsars (iKouveluotou et al. Ill998l . 119991) . The latter class of neutron star s 
with very intense magnetic fields is known as magnetars (IThompson fc Duncan lll993Ul996l ). 
On the other hand, the interior magnetic field could be much higher than t he surface field. 
The l imiting interior field might be estimated from the scalar virial theorem (jLai fc Shapiro 



19911 ). For a typical neutron star mass I.5M0 and radius 15 km, the interior field could be 
as high as ~ 10^^ G. 

Such strong magnetic fields qu antise the motion of charged particles perpendicular to 



the field (ILandau fc Liftshitz 1119771 ). The effects of the phase space modification due to the 
Landau quantisation were studied on the composition and equation of state (EoS) in neutron 
stars extensively. La i and Shapiro extended the Baym, Pethick and Sutherland (BPS) model 



(IBaym et al. Ill971bl ) to the magnetic field case and obtaine d equilibrium nuclei a nd the EoS 



in the outer crust in the presence of strong magnetic fields (iLai &: Shapiro Ill99ll ). The com- 
position and EoS in the core of neutron stars in the presence of strongly quantising magnetic 
fields were investigated w i thin a relativistic field theor etical model by Chakrabarty et al. 



(IChakrabarty et al. 1 119971 : iBandyopadhyay et al. 1 119971 ) . The transport properties such as 



thermal and electr ical conductivities of neutron star crusts in rn agnetic fields were studied 



by several groups (lYakovlev fc Kaminker Ill994l : iHernquist Ill984i ). Recently the magnetised 



(Bavm et al. 


1971a 


) and Harrison- Wheeler EoS for nuclear matter ( 


Nae & Chakrabartv 


2010 


)• 







In the outer crust of a neutron star, neutrons and protons are bound inside nuclei and 
immersed in a uniform background of relativistic electron gas. As the density increases, 
nuclei become more and more neutron rich. Neutrons start to drip out of nuclei at a density 
~ 4 X 10^^ g/cm'^. This is the beginning of the inner crust. The matter in the inner 
crust is made of nuclei embedded in a neutron gas along with the uniform electron gas. 
Further the matter is in /3-equilibrium and maintains charge neutrality. Nuclei are also in 
mechanical equilibrium with the neutron gas. The properties of nuclei in the inner crusts of 
neutron stars in zero magnetic field were studied by different groups. The early studies of 
the inner crust matter were bas ed on the extrapolations of the semiemp irical mass formula 



to the free neutron gas regime (iLanger et al. I Il969l : iBethe et al. I Il970l ) . Baym, Bethe and 



Pethick considered t he reduction of the n uclear surface energy due to the free neutron gas 
in their calculation (IBaym et al. Ill971al ). The study of nuclei in the neutron star crust 
was carried out using the e nergy density of a many body system by Negele and Vautherin 
(INegele fc Vautherin 1 119731 ). With increasing dens i ty in the inner cr ust, unusual nuclear 
shapes might appear there (iRavenhall et al. I Il983l : lOyamatsu 1 119931 ). The properties of 
nuclei in the inner crust were also investigated using a relativistic field theoretical model 
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There are two important aspects of the problem when nuclei are immersed in a neutron 
gas. On the one hand we have to deal with the coexistence of two phases of nuclear matter 
- denser phase inside a nucleus and low density phase outside it, in a thermodynamical con- 
sistent manner. On the other hand, the determination of the surface energy of the interface 
between two phases with good accuracy is needed. It was shown that this problem could 



be solved using the subtr action procedure of Bonche, Levit and Vautherin (IBonche et al. 



1984 . Il985t ISuraud 1119871 ). The properties of a nucleus are isolated from nucleus plus neu 
tron gas in a temperature dependent Hartree-Fock theory using the subtraction procedure. 



This same method was extended to isolated nuclei embedded in a neutron g as (jPe et al 



200ll ) as well as nuclei in the inner crust at zero temperature ( ISil et al. II2002I ). This shows 



that it would be worth studying the properties of nuclei in the inner crust in the presence of 
strongly quantizing magnetic field relevant to magnetars using the subtraction procedure. 

Recently the stability of nuclei embedded in a n electron gas was inve stigated within a 
relativistic mean field model in zero magnetic field ( iBiirvenich et al. 1120071 ). It was observed 
in their calculation that nuclei became more stable against a decay and spontaneous fission 
with increasing electron number density. It is worth mentioning here that the electron num- 
ber density is enhanced in the presence of strong magnetic fields due to Landau quantisation 
compared with the zero field case. The question is whether the nuclear system in the inner 
crust of magnetars would be more stable than those of the field free case. This is the focus 
of our calculation in this article. 

The paper is organised in the following way. In section 2, the formalism for the calcu- 
lation of nuclei of the inner crust immersed in a neutron as well as an electron gas in the 
presence of strongly quantising magnetic fields is described. Results of our calculation are 
discussed in section 3. Section 4 contains the summary and conclusions. 



2. Formalism 

We investigate the properties of nuclei and their stability in the inner crust in the 
presence of strong magnetic fields using the Thomas-Fermi (TF) model. In this case nuclei 
are immersed in a nucleonic gas as well as a uniform background of electrons and may be 
arranged in a lattice. Each lattice volume is replaced by a spherical cell with a nucleus at its 
center in the Wigner-Seitz (WS) approximation. Each cell is taken to be charge neutral such 
that the number of electrons is equal to the number of protons in it. The Coulomb interaction 
between cells is neglected. Electrons are assumed to be uniformly distributed within a cell. 
The system maintains the /3-equilibrium. We assume that the system is placed in a uniform 
magnetic field. Electrons are affected by strongly quantizing magnetic fields. Protons in 
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the cell are affected by magnetic fields only through the charge neutrality condition. The 
int eraction of nuclear ma gnetic moment with the field is negligible in a magnetic field ~ 10^^ 
G Jlattimer et al. Iboooh . 



The calculation below is performed in a zero temperature TF model. In the WS cell, a 
nucleus is located at the centre and immersed in a low density neutron gas whereas protons 
are trapped in the nucleus. However, the spherical cell does not define a nucleus. The 
nucleus is realis ed after subtraction o f the gas part from the cell as shown by Bonche, Levit 
and Vautherin (IBonche et al. I Il984i Il985l ) . In an earlier calc ulation, it was demonstrated 
that the TF formalism at finite temperature gave two solutions ( ISuraud 1119871 ). One solution 
corresponds to the nucleus plus neutron gas and the second one represents only the neutron 
gas. The density profiles of the nucleus plus neutron gas as well as that of the neutron gas 
are obtained self-consistently in the TF formalism. Finally the nucleus is obtained as the 
difference of two solutions. This formalism is adopted in our calculation at zero temperature 
as described below. 

The nucleus plus gas solution coincides with the gas solution at lar ge distance leading to 
the d efinition of the thermodynamic potential (^at) of the nucleus as (IBonche et al. Ill984j . 
1985h 

= ^NG — 5 (1) 

where Qng is the thermodynamic potential of the nucleus plus gas phase and is that of 
the gas only. The thermodynamic potential is defined as 



(2) 



q=n,p 



where Hq and Aq are the chemical potential and number of q-th species, respectively. The 
free energy is given by 

J^iuh, Yp) = J['H + e, + Eeldr , (3) 

where "H is nuclear energy density functional, Ec is the Coulomb energy density and Se is 
the energy density of electrons. The free energy is a function of average baryon density 
{rih) and proton fraction (Yp). The nucle ar energy densi t y is calculated using the SKM 



nucle on-nucleon interaction and given by (ISil et al. 1 120021 : iBrack et al. I Il985l : IStone et al. 
2003h 
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The first two terms of the nuclear energy density are kinetic energy densities of neutrons 
and protons, respectively. The third term originates from the zero range part of the Skyrme 
interaction whereas the last term is the contribution of the density dependent part of the 
nucleon-nucleon interaction. The effective mass of nucleons is given by 
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where total baryon density is p = p„ + Pp. 

The Coulomb energy densities for the NG and G phases are: 



r — r 

^2 



r — r 



— (ra^(r') - ne)dr' 



r — r 



— (n^(r') - ne)dr' 



(6) 



where nFj^Q and Uq are proton densities in the nucleus plus gas and only gas. Here the 
coulomb energy densities (Ec) represent the direct part. We do not consider the exchange 
part because its contribution is small. 

So far the formalism described above is applicable for the zero magnetic field. However, 
we study the effects of magnetic fields on electrons which , in turn, influence the properties 
of nuclei in the inner crust. The Coulomb energy density and the energy density of electrons 
in Eq. (3) would be modified in strongly quantising magnetic fields. In the presence of a 
magnetic field, the motion of electrons is quantized in the plane perpendicular to the field. 
Protons in nuclei would be influenced by a magnetic field through the charge neutrality 
condition. We take the magnetic field {^) along Z-direction and assume that it is uniform 
throughout the inner crust. If the field strength exceeds a critical value Be = ml/e ~ 
4.414 X 10^'^G, then electrons become relativistic (iLai I l200l[ ). The energy eigenvalue of 
relativistic electrons in quantizing magnetic field is given by 



2z/R 



1/2 



(7) 
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where pz is the Z-component of momentum, B^, = B/Bc, v is the Landau quantum number. 
The Fermi momentum of electrons, p/^, is obtained from 

where V"'^(r) is the direct part of the single particle Coulomb potential. 
The number density of electrons in a magnetic field is calculated as 

eB " 



27r^ 



(9) 



u=0 



Here the spin degeneracy is = 1 for the lowest Landau level (z/ = 0) and = 2 for all 
other levels. 



The maximum Landau quantum number {vmax) is given by 

_ {yi, + V\r)f -ml 
2eB 

The energy density of electrons is obtained from. 



(10) 











We minimise the thermodynamic potential in the TF approximation with the condition 
of number conservation of each species. The density profiles of neutrons and protons with 
or without magnetic fields are obtained from 
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(Sil et al. 


2002; 


De et al. 


2001) 



(37r2)i-^(nV) 



'+Vng + V^nin 



NGK"'NG^ 



2 2 



2m* 



(12) 



(13) 



where m* is the effective mass of q-th species , V^q and V p are the single par ticle potentials 



of nucleons in the nucleus plus gas as well as gas phases (IBrack et al. Ill985l ). On the other 
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hand, V^q and Vq are direct parts of the single particle Coulomb potential corresponding 
to the nucleus plus gas and only gas solutions and both are given by 

V^r) = J KG(rO - ne] ^^^^dv' . (14) 

The average chemical potential for q-th nucleon is 

^'^tJ + V^oir) + V^Gir)]p%Gir)dr (15) 

where Aq refers to Nceii or Zceii of the cell which is defined by the average baryon density rib 
and proton fraction Yp. The /9-equilibrium condition is written as 

l^n ^ fJ.p + IJ,e ■ (16) 
The average electron chemical potential in magnetic fields is given by 

^ie^[pM' + ml{l + 2uB^)Y/'-<V'^ir) > , (17) 

where < V^lr) > denotes the average of the single particle Coulomb potential. 
Density profiles of neutrons and protons in the cell are constrained as 

Nceii = j Pn''{r)dv , (18) 

where Nceii and Zceii are neutron and proton numbers in the cell, respectively. 

Finally, number of neutrons (iV) and protons (Z) in a nucleus with mass number A — 
N + Z are obtained using the subtraction procedure as 

N = J [P^''{r)-P^{r)]dr. (19) 

3. Results and Discussion 

We find out the equilibrium nucleus at each density point minimising the free energy of 
the system within a WS cell maintaining charge neutrality and ^-equilibrium. The variables 
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of this problem are the average baryon density (n^) , the proton fraction (Yp) and the radius 
of a cell {Rc). For a fixed value of n^, Yp and Rc, the total number of nucleons (Aceii) is given 
by Aceii = VceiiUh where the volume of a cell is Vceii = 4/37ri?^. The proton number in the 
cell is Zceii = YpUhVceii and the neutron number is Nceii = Aceii — Zceii- Now we obtain density 
profiles of neutrons and protons in the cell using Eqs. (13) and (18) at a given average baryon 
density and Yp. Consequently, we calculate chemical potentials of neutrons and protons and 
free energy per nucleon. Next we vary the proton fraction, calculate chemical potentials and 
density profiles and obtain the ^-equilibrium in the cell. Finally we adust the cell size (Rc) 
and repeat the above mentioned steps to get the minimum of the free energy. These values 
of Yp and Rc are then used to calculate neutron and proton numbers in a nucleus at an 
average baryon density corresponding to the free energy minimum with the help of Eq. (19). 
This procedure is repeated for each average baryon density. 

The minimum of the free energy originates from the interplay of different contributions. 
The free energy per nucleon is given by 

F/A = ejv + eiat + e^ie ■ (20) 

The nuclear energy including the Coulomb interaction among protons is denoted by e^v, eiat is 
the lattice energy which involves the Coulomb interaction between electrons and protons and 
the electron kinetic energy is Ceie- The free energy per nucleon in the presence of magnetic 
field S = 4.414xlO^^Gis shown as a function of the cell size for an average baryon density 
rib = 0.008 fm ~^ in Fig. 1. We note that the nuclear energy increases with Rq. On the other 
hand, the lattice energy and electron kinetic energy both decrease with increasing cell size. 
The competition of e^v with the sum of eiat and Ceie determines the free energy minimum. 
The cell radius corresponding to the free energy minimum is 32.1 fm for the zero field case 
(not shown in the figure) and 31.9 fm for B = 4.414 x 10^^ G. The corresponding proton 
fraction for B = AAU x 10^^ G is 0.03. 

In Figure 2, the cell size corresponding to the free energy minimum is plotted as a 
function of average baryon density for magnetic fields B = 0, 4.414 x 10^*', 10^'' and 4.414 x 
10^'^ G. For magnetic fields B < 10^^, several Landau levels are populated by electrons. 
Consequently we do not find any change in the cell size in the magnetic fields compared 
with the zero field case. However we find some change in the cell size for S = 10^^ G when 
only the zeroth Landau level is populated by electrons for ni, < 0.004, fm~^ whereas first 
two levels are populated in the density range 0.005 to 0.015 fm"^. However the cell size is 
increased due to the population of only zeroth Landau level in the presence of the magnetic 
field B = 4.414 x 10^''G compared with the zero field case. The size of the cell always 
decreases with increasing average baryon density. 

The proton fraction in the presence of magnetic fields is shown as a function of average 
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baryon density in Fig. 3. Protons in nuclei are affected by the Landau quantisation of 
electrons through the charge neutrality condition in a cell. For magnetic fields B < 10^^ G, 
the proton fraction is the same as that of the zero field case over the whole density range 
considered here. We find some changes in the proton fraction below = 0.015 fm~^ when 
the field is 10^^ G. Though electrons populate the zeroth Landau level for nb < 0.004 fm~^, 
the proton fraction decreases below the corresponding proton fraction of the zero field case. 
However, for magnetic field B = 4.414 x 10^'' G, the proton fraction is strongly enhanced 
due to the population of the zeroth Landau level for < 0.04 fm~^ compared with the zero 
field case. 

The density profile of neutrons in the nucleus plus gas (NG) as well as gas (G) phases 
corresponding to rih = 0.02 fm~'^ with and without magnetic fields are exhibited as a function 
of distance (r) within the cell in Fig. 4. The solid line denotes the zero field case where as the 
dashed line represents the density profile with the field B = 4.414 x 10^^ G. The horizontal 
lines imply the uniform gas phases in both cases. The proton fraction is 0.040 for the magnetic 
field case whereas it is 0.022 for the zero field case. Though neutrons are not directly affected 
by the magnetic field, the neutron chemical potential is modified through the /3-equilibrium 
due to Landau quantisation of electrons. Consequently, the number density of neutrons is 
altered. We find that the neutron density is higher in the gas phase for the zero field case 
than that of the situation with the magnetic field. We show the subtracted density profiles 
of neutrons with magnetic field B = 4.414 x 10^^ G in Fig. 5. The neutron density profiles 
in the nucleus phase with and without magnetic field are different. Further we note that 
less number of neutrons drip out of a nucleus in the presence of the magnetic field than the 
situation without the field. This may be attributed to the shift in the /3-equilibrium in strong 
magnetic fields. We find a similar situation in the calculation of t he outer crust in magnetic 



fields that the neutron drip point is shifted to higher densities (INandi fc Bandyopadhyay 



20101). 



Now we know the density profiles of neutron and protons in the nucleus plus gas phase 
as well as in the nucleus at each average baryon density. We immediately calculate the total 
number of neutrons and protons in the nucleus plus gas phase and in a nucleus using Eqs. 
(18) and (19). We show total number [Aceii) in a cell for magnetic fields B = 4.414 x 10^®, 
10^^ and 4.414 x 10^^ G with average baryon density in Fig. 6. The dotted line denotes 
the zero field case. In all cases, Aceii growing with the density reaches a maximum and 
then decreases. Such a trend w as observed in the calculat ion of Negele and Vautherin in 



the absence of a magnetic field (INegele fc Vautherin Ill973l ). We note that our predictions 
for B = 4.414 x 10^^ G do not change from the field free results because a large number of 
Landau levels is populated in that magnetic field. For mag netic field B = 10^^ G, the total 
number of nucleons decreases in the density regime 0.005 - 0.02 fm~^ compared with the 
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corresponding results of the field free case. This can be understood from the behaviour of 
the cell size around that density regime in Fig. 2. For B = 4AU x 10^'^ G, the zeroth Landau 
level is populated by electrons for densities < 0.04 fm^'^ This modifies the /3-equilibrium 
and the charge neutrality conditions which, in turn, impact the size of the cell and the total 
number of nucleons in a cell. This effect is pronounced in the case of -B = 4.414 x 10^^ G. In 
this case, Aceii is significantly reduced compared with the zero field case for densities < 0.04 
fm~^. 

We obtain neutron (A^), proton (Z) and total nucleon numbers (A) in the nucleus at 
each average baryon density following the subtraction procedure. Total nucleon and proton 
numbers are shown in Fig. 7 for the above mentioned magnetic fields. When the magnetic 
field is 4.414 X 10^^ G or more. It is noted that our results in certain cases start oscillating from 
the field free results. This may be attributed to the fact that the population of Landau levels 
jumps from a few levels to zero in the above mentioned fields as baryon density decreases from 
higher to lower values. In contrast to Fig. 6, we find total nucleon and proton numbers inside 
the nucleus at each density point beyond 0.002 upto 0.04 fm~^ are significantly enhanced in 
case of -B = 4.414 x 10^'' G compared with the field free case as well as other magnetic fields 
considered here. This clearly demonstrates that more neutrons are inside the nucleus in the 
presence of strong magnetic fields > 10^^ G than in the gas phase in that density regime. 
This is opposite to the situation in the zero magnetic field. This can be easily understood 
from the density profiles with and without magnetic fields in Fig. 4 and Fig. 5. 

We plot the free energy per nucleon of the system with average baryon density in Fig. 8. 
Our results for B = 4.414 x 10^^ G do not change much form the field free results. However, 
for B — 10^^ G, the free energy per nucleon is reduced at lower densities (< 0.004 fm~^) as 
compared with the field free case. We find more pronounced reduction in the free energy 
per nucleon in the field B — 4.414 x 10^^ G almost over the whole density regime considered 
here. 



4. Summary and Conclusions 

We have investigated properties of nuclei in the inner crust and their stability in the 
presence of strong magnetic fields ~ 10^^ or more. Nuclei are immersed in a neutron gas and 
uniform background of electrons. We have adopted the SKM* interaction for the nuclear 
energy density functional and studied this problem in the Thomas-Fermi model. Electrons 
are affected through Landau quantisation in strong magnetic fields because much less Landau 
levels can be occupied in these cases. Consequently, electron number density and energy 
density are modified in strongly quantising magnetic field and the /3-equilibrium condition 
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is altered compared with the field free case. The enhancement of electron number density 
in magnetic fields > 10^^ G due to the population of the zeroth Landau level leads to 
enhancement in proton fraction through the charge neutrality condition. We minimise the 
free energy of the system within a Wigner-Seitz cell to obtain the nucleus at each average 
baryon density. In this connection we used the subtraction procedure to obtain the density 
profiles of a nucleus from the nucleus plus gas and only gas solutions at each average baryon 
density point. We note that less number of neutrons drip out of a nucleus in the presence 
of strong fields than the situation without magnetic field. This results in larger mass and 
proton numbers in a nucleus in the presence of magnetic field > 10^'' G compared with the 
corresponding nucleus in the field free case. Further the free energy per nucleon of the system 
is reduced in magnetic fields > 10^'' G. 

Magnet ars might eject crustal matter due to tremendous magnetic stress on the crust 
( iGelfand et al. 1120051 ). The ejected matter of the inner crust might expand to much lower 
densities. The decor npressed crustal matter has l ong been considered as an important site 
for r-process nuclei (ILattimer fc Schramm 1119771 : lArnould et al. 1120071 ). It would be worth 
studying the r-process in the decompressed crustal matter of magnetars using the results of 
our calculation as an input. 



RN and DB thank the Alexander von Humboldt Foundation for the support under the 
Research Group Linkage programme. 
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Fig. 1. — Free energy per nucleon is plotted with cell size for average baryon density rit, — 
0.008 fm-3 and magnetic field B = 4AU x 10^^ G. 
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Fig. 2. — Cell size corresponding to the free energy minimum is shown as a function of 
average baryon density for magnetic fields B = 0, 4.414 x 10^^, 10^'' and 4.414 x 10^'' G. 
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Fig. 3. — Proton fraction is plotted with average baryon density for magnetic fields B — 0, 
4.414 X 10^^ 10^^ and 4.414 x 10^^ G. 
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Fig. 4. — Neutron density profile in the nucleus plus gas phase with magnetic field B — 
4.414 X 10^^ G (dotted fine) and without magnetic field (sohd fine) at an average baryon 
density rib — 0.02 fm~^. Horizontal fines denote gas phases in both cases. 
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Fig. 5. — Neutron density profile in the nucleus plus gas phase and after the subtraction of 
the gas phase with magnetic field B — 4.414 x 10^^ G for same values of average baryon 
density and proton fraction as in Fig. 4. 
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Fig. 6. — Total number of neutrons and protons in a cell (Aceii) is plotted as a function of 
average baryon density for magnetic fields B = 0, 4.414 x lO^'', 10^^ and 4.414 x 10^^ G. 
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Fig. 7. — Total nucleoli and proton numbers in a nucleus are plotted as a function of average 
baryon density for magnetic fields B = 0, 4.414 x 10^*^, 10^'^ and 4.414 x 10^^ G. 
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Fig. 8. — Minimum Free energy per nucleon of the system is shown as a function of average 
baryon density for magnetic fields fields B = 0, 4.414 x 10^*^, 10^^ and 4.414 x 10^'^ G. 



